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We present an analysis of the SU(3) symmetric model of the strongly interacting three component 
Fermi gas in the continuum space using quantum Monte Carlo techniques. Three body effects pre- 
dominate in the regime of interaction strength beyond that of threshold of the three particle bound 
state. However, we find that there is an interval of the interaction strength where the SU(2)®U(1) 
broken symmetry superfluidity is possible. For a strong enough interaction, the SU(3) symmetry 
is restored and the superfluidity is suppressed. Within the interval of the broken symmetry, we 
also find that on average the particle pairs belonging to the species with superfluid pairing remain 
further separated than those without the superfluid pairing correlation. 

PACS numbers: 03.75.Ss, 05.30.Fk, 21.65.-f, 31.15.A-, 31.15.xt 



I. INTRODUCTION 



A pair of fermions occupying different internal states can interact in s-wave scattering. Fermi atoms can be loaded 
into different internal spin projection states to produce interacting degenerate Fermi gas with T m 0. Comprehensive 
reviews of the recent theoretical and experimental advances in the study of the dilute Fermi gases can be found in 
the Ref. [l], . The fermionic alkali atoms can occupy hyperfine states that result from the coupling of the nuclear 
angular moment / with the electronic spin S. In the presence of a strong magnetic field, the electronic degree of 
freedom is polarized and the relevant internal degrees of freedom are determined by the state of the nuclear angular 
momentum. In the case of 6 Li, the nuclear angular momentum I — 1 allows three different projections. By combining 
magnetic as well as optical trapping fields, constraints on the trappable hyperfine states are largely removed. 

For a Fermi gas at low density, an expansion of the ground state energy in terms of a small parameter ak F is known 
[USUI- In general, when we have a gas of fermions occupying s different states with the same partial densities, we 
call it the degeneracy s Fermi gas. The ground state energy has an expansion in powers of ak F 
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is the Fermi energy, k F 



where we identify the zeroth order term as the free Fermi gas energy E F q = % e F- £f — . 
the Fermi momentum and a the s-wave scattering length. E F q is independent of s and it is used as the unit of energy 
throughout this article. The terms that depend on the interaction potential range R are to appear at higher orders 
0- For s = 2 Fermi gas, the dependence on R can be eliminated by taking the limit R/ro — > 0, where tq is the 
average inter-particle distance (~7ITq/0 — 1 with p = density). The intermediate regime where R « r ~ j- << \a\ 
is of particular interest for s = 2 Fermi gas. This is also called the unitarity regime. In this regime, the mean free 
path of the atoms A = becomes much shorter than ro as the cross section a diverges. Examples of the systems in 

this regime are found in the dilute gases of 6 Li and 40 K atoms at Feshbach resonances close to the zero temperature. 
The magnitude of the s-wave scattering length \a\ can be ~ 1000A at the Feshbach resonance while the interaction 
range of the van der Waals forces R ~ 10 — 100A Another example is that of the neutron gas. The neutron-neutron 
interactions by strong force have a ~ — 18.8/m while R ~ lfm. At a — > ±oo, zero energy two particle bound state 
appears. However, it is known that the ground state of the many body system has a positive energy per particle 
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E /{NE FG ) w 0.40 - 0.44 @, [§]. Ground state properties for s = 2 Fermi gas were studied@, 0, using the 
Quantum Monte Carlo for the unitarity as well as other regimes of interaction. The generalization of the low density 
expansion (Eq. [IJ to s > 3 becomes troublesome as the limit of R — > cannot be taken. This is expected since for 
the three particle systems the Efimov effect predicts an effective three particle attractive interaction Q of the form 

s 2 n 2 

~ — 2 mR: i wnere so is a universal constant. Thus, the minimum set of parameters to describe the s = 3 system consist 
of akp and Rkp. In general, the mean field treatments of the three component Fermi gas [IH [IH [3 do not account 
correctly the three particle physics. 

Recently, a stable three component degenerate Fermi gas [l5| has been created experimentally. Broad and close 
lying Fcshbach resonances make the strong and attractive interactions among the atoms in the different internal states 
possible [ll| . For simplicity, we label the atoms occupying these three different states by a coloring scheme: Green (G), 
Red (R) and Blue (B). The channel dependent s-wave scattering lengths (clgr, &GB, and club) have been measured 
experimentally fl7j. It is experimentally difficult to achieve simultaneously strong interactions in all of the channels. 
However, in the present work we assume a simplified SU(3) symmetric model where a — agjj — o-gb — o-rb- Also, the 
mass is assumed to be the same for all the components. We leave non-SU(3) symmetric cases to future study. Three 
component Fermi gas problem is also relevant in relation to the color superconductivity of the quark matter [l8| . 

In this article, ab initio Monte Carlo results of the three component Fermi gas are presented. The fermions 
interact pairwise, by a short but finite range attractive interaction potential. The possibility of the SU(2)®U(1) 
broken symmetry ground state [H, [T|| is considered. Here, only two components participate in the superfluid pairing 
while the third component remains in the normal phase. At weak interaction strengths, the superfluid pairing is 
exponentially suppressed with the quasiparticle gap A/Tp ~ e 7T ^ 2akF \ However, we find that there is an interval 
of the interaction strength where the broken symmetry pairing has noticeable effects in the ground state energy, the 
quasiparticle gap, the pair distribution functions, etc. When the strength of the interaction is further increased, the 
SU(3) symmetry is restored by the predominant three body effects. We present our analysis using the dimensional 
arguments for the three particle system (Sec. |Hj and the quantum Monte Carlo method for the three component 
Fermi gas (Sec. 111111 Vp . 



II. THREE PARTICLE BOUND STATE AND SCALING BEHAVIOR 



L.H. Thomas [2(| noted that in the nucleus of tritium ( 3 H) which has two neutrons and one proton, the binding 
energy has no lower bound if we assumed finite negative s-wave scattering length for the proton-neutron interaction 
and took the interaction range to zero. Thus, three particle binding energies depend on the range of the potential 
and the ground state energy diverges in the limit of zero interaction range. In addition, as the pairwise interaction 
approaches the resonance (a — > — oo), infinite number of shallow three particle bound states (called trimers or trions) 
appear one after another. They are known as the Efimov states [2l|, [H, HH, H3 ■ This property is strikingly different 
from the two particle case. Consequently, qualitatively different behavior of the three component gas is expected to 
emerge in comparison with the two component Fermi gas. For our discussion, we consider a generic three particle 
Hamiltonian where the particles interact pairwise 

= E V '+^ E Vk(r«). (2) 

i£{G,R : B} i<je{G,R,B} 

There are two positive parameters in the Hamiltonian: strength vq and range R of the interaction. We assume a 
negative dimensionless function Vr(t) < while vq > has dimension ~ length" 2 . Vr(t") is solely parametrized by 
R. The s-wave scattering length a can be used equivalently instead of the strength vq. 

In the two particle systems the bound state threshold is at the resonance (a^ = ±oo) independent of the finite 
R. Given a value of i?, it is easy to see that three particle threshold a§ can be such that a§ < but not equal to 
—oo. This can be done using an ansatz as the trial wave function; for example, a function of the form ^s^riali^-) = 
f{rGR)f{'i'GB)f{'i'RB) with the variational f(r). In this case, it suffices to provide a trimer state upper bound in 
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energy to prove the existence of a trimer state with a| in the interval (— oo, 0) for a given R. The fact that is 
negative and ^ —00 is the starting point for the analysis on the scaling behavior of the length unit. When we 
consider rescaling of the length by taking R — > aR and a —> aa (we will always assume < a < 1 from now and on), 
consistent scaling behavior with (H3) — > (Ti.3}/a 2 is expected (we will justify in the next paragraph). Here, it can be 
seen easily from the zero energy scattering solution that the s-wave scattering length a scales analogous to another 
length quantity R. By making a — > + we expect R — > + , a§ — > 0~ and (H3) — » —00 for any a that belongs to the 
interval (-co, a§ — > _ ). This means that, for a zero range interaction, the trimer state is possible for any attractive 
pairwise interaction of nonzero strength. 



Now, we consider formally the scaling behavior of the energy for the system of three particles in vacuum. We 
assume a value of R such that trimer state is allowed with a strength Vq corresponding to a§ < 0. The three particle 
Schrodinger equation with the usual notation X = {x<3,x fl ,x B } and T%j |x^ JCj J is 



W 3 *(X) = EV(X) 



(3) 
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and tf s (X s ) = tf(X). Then Vf* s (X s ) = ^V, 2 *(X) and 



Vu(rij)^> (X.) . Thus, after scaling all the length quantities by a, the Eq. [3] becomes 
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As a result, after length scaling \1/ S (X S ) is the solution with the eigenvalue E s = We arrived at this property by 
using the dimensional arguments alone. As good illustrative examples of this scaling behavior, the following cases are 
mentioned: When a G (—00, a|) for a given R, there is trimer state with energy E < 0. Then, we can scale length by 
an overall factor a but keep a s constant (that is, increase \a s \ to match \a\). It is obvious that taking a — ► + causes 
E s (a) to collapse rapidly to — 00. This means that given a trimer state in vacuum, when R — > at a fixed non-zero 
value of a (< 0), E has to go to —00. This is in agreement with the above mentioned Ref. [2(|. As another example, 
if we had initially a £ (a§ , 0) , the interaction potential is not strong enough to allow trimer state and the total energy 
of the particles E = in the vacuum. In this case, length scaling leaves the particles unbound with E s (a s ) = 0. Thus, 
unbound particles remain unbound even after length scaling and the collapse does not occur. 



For the comparison purpose, let's also consider the scaling behavior of a pair of particles. We can see that the 
scaling of length does not produce the collapse as in the three particle case. For the s-wave scattering length a such 
that l/a < 0, there is no bound state and the energy of the pair remains zero {E = 0) in the vacuum at any length 
scale. Then, we consider l/a > regime (usually called BEC regime). Using the same scaling analysis, we also have 

E s (a s ) — • In this case, we can solve exactly the contact interaction (R — 0) problem by replacing the potential 
by the boundary condition = — i where u(r) is the radial wave function of the pair. The solution for the radial 
wave function is u(r) ~ e~~ with E pa i r {a) — —-^i- This energy is finite unless a — > + . According to this solution, 
the scaling behavior of energy E p°-"-(° ia ) — j s ^ ne same as the result we obtained from the dimensional analysis, 

Impair ) OL 

Es ^°^ = 7^2. Thus, for a pair the energy scaling relation becomes identical to the bound state energy in the strong 
coupling limit l/a s —> +00. This is qualitatively different from the three body collapse at the weak coupling limit 
X/a s — > — 00 because of the simultaneous scaling of R to zero. 
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III. MANY PARTICLE GROUND STATE 



For the study of many body systems at finite density, we use ab initio stochastic method known as Fixed Node 
Green's Function Monte Carlo (FN GFMC). In general, we take a trial wave function <Jv that obeys antisymmetry 
upon the exchange of identical fermions and let it evolve in the imaginary time restricted to a definite sign domain 
given by the nodal surface of the trial wave function itself. If the nodal structure is correct, we get the exact ground 
state. Otherwise, we get an approximate ground state and an energy upper-bound. The implementation of this 
method is explained in detail elsewhere [1, [ll| . 



We consider a system with 6^8 particles of each color (Green, Red, Blue), so that 18 < Ntotai < 24. The particles 
are contained in a finite box with the periodic boundary conditions at the walls to simulate the uniform matter. For 
many particle systems, we cannot rescale the length without changing the density. The scaling behavior analyzed 
in the previous section is only applicable to the few particle systems in vacuum. In the degeneracy two case, one 
dimcnsionless product akp uniquely determines the system. The parameter R can be pushed in principle to the 
zero limit and eliminated from the description of the system. In practice, small but finite values of R/r Q << 1 were 
assumed [ll[ as long as the results were converged within the statistical errors to the RRf — > limit. However, 
from the scaling behavior analysis of the previous section, it becomes clear that for the degeneracy three Fermi gas we 
need to keep RkF finite in order to avoid local trimer instability. Thus, we need both akF and RkF to fully describe 
the three component Fermi gas. We keep Rkp — 0.32 in the s — 3 case which is the same value as in the s = 2 case 
(i?/ro ~ 0.1, thus converged to the R/ro — > limit). However, we should keep in mind that for s — 3, this particular 
value of R is not the limit of R — > 0. We consider the cases with small deviations from the balanced partial densities. 
We explore the possibility of the superfluid ground state within a certain interval of the interaction strength. We also 
analyze the qualitative behavior of the pair distribution functions (see Sec lIV[) . The many body SU(3) symmetric 
Hamiltonian is 



H = E v * + v ° E vh(r«)(i - 



(5) 



where Cj is the color of i-th particle. Only pairs of different color particles interact. Vr[t) is the dimensionless core 
of the Poschl- Teller potential: Vr[t) = — cosh -i Hf/m an< ^ v ° ~ ^W- v ° can ^ e a djusted to get the desired scattering 
lengths. We can impose different nodal restrictions to the solution by using different trial wave functions. We can 
estimate the energy by using SU(3) symmetric Slater trial wave function 
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(6) 



where the factors represent the normal states (given by Slater determinants of the plane wave orbitals) of different 



color species. It was also shown[12| that the pairing fields (A Gfl , A GB , A# B ) with A Qj a ~ Sk( c k,. 
mapped into (Aq, 0, 0) with the constraint ^ Q a A^ = Aq. This is analogous to the analysis of Rcf. 
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and Leggett where the ground state should allow one normal phase component. Thus, we consider the SU(2)® U(l) 
broken symmetry pairing ground state with BCS pairing for two of the Fermi components while the third component 
remains in the normal phase. The corresponding nodal structure is given by the trial wave function 
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Here, we assume that blue species remains in the normal phase (represented by ^ fg,b), while between the green and 
red species there is pairing correlation (represented by ^ bcs,Gr)- In the last line of Eq. [71 we assumed fixed number 
projection of the green and red particles. 



The complete trial wave function with Jastrow-like factor can be written as 

*v,fg - I| [fs{r i0 )8 Cl , C] + f D (nj)(l - S Ct . Cj )] * FG 



(8) 



l<3 
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for the SU(3) symmetric Slater wave function (Eq. [6]). Here, we considered separately the correlation that exists 
between the same color particles /g(r) (= /ggM = fRR.{ r ) — /bbW) from the correlation that exists between 
different color particles /o(r) (= fcR{ r ) = /gbM — /rb(V)). Usually fs(r) is qualitatively different from /o(r). 
fs(r) is analogous to /tt( r ) °f f ne degeneracy two case and includes Pauli exclusion principle (/s(0) = 0). The 
particular shapes of the correlation functions do not change the GFMC energies. However, they are optimized in 
order to minimize statistical errors and to have the converged pair distribution functions g(r) « gtriaifj 1 ) ss gGFMcir) 
(see Ref. [25(). We can see that the optimized foiTij) deviates largely (more peaked at r ~ 0) from the one obtained 
by using the LOCV equations [8(. This is due to the strong three body effects even when the pairwise interactions 
are relatively weak. 

Analogously, the complete trial wave function with the SU(2)(g>U(l) broken symmetry pairing correlation (Eq. [7]) 

is 

^VM-BCS = Y\[fs{nj)8c i ,c ] +fGR.(r lj )5 CuG 5 Cj jt + ... 

fGB{r i] )(5 Cz ,G + Sc^^Sc^B^bs-BCS ■ (9) 

This wave function can give better pair distribution functions as the optimization of the correlation functions can 
be carried out separately for /crM and /gb(^) = /abM- Then, we can demonstrate that Jgr(0) < /gb(0) as 
well as <?Gi?(0) < #gb(0) (see the discussion of the Fig. [5] in the Sec. IIV|) . The optimum nodal structure is tried as 
ui = {0.2, 0.1, 0.02, O.Of , 0.01} with non zero short range function /3(r) parametrized by b — 0.44 (see Ref. [1, [ll| for 
the definitions of these parameters). These parameters are identical to those of the s — 2 case when > (called 

molecular or BEC regime) . This trial wave function assumes arbitrarily one of the three possible broken symmetry 
pairing states 

|1) = ^FCB^BCStGR 
|2) = ^FG^^BCS.GB 

|3) = ^FG,G^BCS,RB- (10) 

These states are degenerate in energy and the broken symmetry can be assumed by choosing one of these states 
without loss of generality. The GFMC energies using the trial nodes given by the wave functions of Eq. [8] and Eq. [9] 
are summarized in the Fig. [TJ while the pair correlation functions are presented in the Fig. [2l 



IV. DISCUSSION OF THE RESULTS 



The results of three component Fermi gas energies are summarized in the Fig. Q] The pairing correlations produce 
noticeable effects at —1 < 1/akF < —0.5 (compare the triangles with the circles of the same figure). Close to the 
ag, the pairing effects are small. Here, the energy of the SU(2)(g)U(l) broken symmetry state is not distinguishable 
within the error bars from that of the SU(3) symmetric state (Slater) . The system is found stable in the regime of 
interaction considered in this work (—1.3 < 1/akp < —0.3 and Rkp = 0.32). 

We also found that it is possible to see the effects of pairing in the g(r)'s. In the plot a) of Fig. [5] (akp — —0.75), 
no difference can be seen in the <7(r)'s because the pairing is weak. However, where the pairing is relevant (in terms of 
the energy) we notice that g G R.(0) < 3gs(0) = g_R B (0) < gsiater(0) (where g S iater(r) is the pair distribution obtained 
with the non-pairing Slater wave function) and the symmetry is broken (see the plots c) and d) of the Fig \2§- Green 
and red particles have superfluid pairing correlations, so one may have naively expected gGfl(0) > <7gb(0) = .9/?,s(0), 
but the opposite is found to be true. According to this, the relative distances satisfy on average: tqb — Trb < r G R- 
The interpretation we can give is that once a Cooper pair is formed, the third particle feels much stronger attraction 
toward the center of the mass of the Cooper pair. In fact, we can approximately estimate that the potential strength 
between the center of mass (CM) of green-red pair and the blue particle is enhanced from vq to |vo, while green and 
red particles interact with the strength vq. This comes from the observation that the zero energy scattering equation 
between the CM of green-red and blue is —^- u "( r ) + VQVn(r)u(r) = with m r = |m. From the plot d) of the 
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FIG. 1: (color online) Comparison of the degeneracy two (s = 2) and three (s = 3) results for RkF = 0.32 in both cases. ELenz 
(Eq. [TJ gives a good estimate of the normal state energy for s = 2 and 1/akp < —1. However for s = 3, the match is poor 
with the GFMC results as the R dependent terms and three particle effects predominate. In the region 1/afci? ^ —0.3, the 
SU(2)®U(1) broken symmetry pairing state is shown to be the ground state for s = 3 Fermi gas. However, for —0.3 < l/afci? 
the SU(3) symmetry is restored(see also Fig. [4]) and the superfluid pairing suppressed. This behavior for s — 3 is quantitatively 
different from the s — 2 case, where with the increasing interaction strength (increasing 1/akp), the superfluid pairing state 
(crosses) is always the favored ground state compared to the state without superfluid correlation (circles). The size of the 
symbols approximately correspond to the error bars. 
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FIG. 2: (color online) Pair distribution functions at afci? = —0.75(a), —1.0(b), —1.3(c) and —2.0(d). We define pair distributions 
9a/3(r) ~ ( y, 8(r — ri a ,ja))- The angled brackets imply taking thermal average of the possible configurations. g a p(r) has 

itj ,3f3 

boundary condition lim g a p(r) = 1. The line corresponds to gcRir), that is, the pair distribution between the green and red 

r — >oo 

species while the dashed line corresponds to the distribution between the green and blue species gGB{r). Both gGii(r) and 
gGB(r) assume SU(2)®U(1) symmetric state (Eq. [9|. The dotted line represents the pair distribution without superfluidity 
gsiater(r) calculated with the SU(3) symmetric Slater wave function (Eq. [8}. The gsiater(r) is independent of the colors as long 
as they are different. All g(r)'s are calculated with the optimized f(r) functions. In the plots c) and d) where the superfluid 
pairing is strong, gGFt(0) < ffGfl(O) = 3hs(0) < gsiater(0)- Thus, on the average the particles that have superfluid pairing 
correlation remain further apart than the particles that have no superfluid pairing correlations. 
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FIG. 3: (color online) Energy excitations with unbalanced particle numbers (allows broken pairs) at akF = —1.3. The circles 
correspond to the set of particle numbers {N G , Nr, Nb} = {6, 6, 6}, {7, 6, 6}, {7, 7, 6}, {7, 7, 7}, {8, 7, 7}, {8, 8, 7}, and {8,8,8}. 
The squares are {6, 6, 7}, {6, 6, 8}, {7, 7, 8}, and {7, 7, 9} cases. The triangles represent the configurations {7, 6, 7} and {8, 7, 8}. 
The size of the symbols approximately correspond to the error bars. 
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FIG. 4: (color online) Imaginary time evolution of the overlap of the wave functions defined as (2|e~ WT |l) = (3|e^ T |l) ~ 
(^/_FG|e _Wr |l) are shown. Each time step is 5t = 3.6 x 10 4 E n p ■ The states |1), |2) and |3) are as defined in Eq. 1101 while 
I^fg) is as defined in Eq. [8] We can see that as the interaction strength is increased, the symmetry is restored with the 
overlaps approaching 1. 



Fig. [2J it is still not conclusive whether the symmetry could be restored (<?gr(0) = 5gb(0) = guB (0)) in the regime 
of strong three particle correlation —0.5 < 1/akp. In that regime, the analysis is hindered on the practical ground: 
it becomes increasingly harder to obtain reliable estimates of g(r) keeping the statistical errors small. However, the 
Fig. [1] and Fig. [4] show clearly the tendency toward the restoration of the SU(3) symmetry entering into a regime 
completely dominated by three particle physics with negligible paring correlations. 



Using the GFMC technique, energy gaps can be calculated by allowing variations in the numbers of the green and red 
particles while keeping the number of blue particles constant. In this way, we allow breaking of the superfluid pairs. 
We can consider sets of {N G ,N R ,N B } = {6, 6, 6}, {7, 6, 6}, {7, 7, 6}, and {7, 7, 7}, {8, 7, 7}, {8, 8, 7}. We observe the 
usual odd-even staggering of the ground state energy (Fig. [3]). Furthermore, we can allow the variations 6 < Nq < 8, 
6 < N R < 8, and 6 < Nb < 8 and calculate the excitation energies. The lowest energy excitations at akp = —1 and 
— 1.3 are found with the momentum k = quasiparticlc. In the case of s = 2, quasiparticles with zero momentum 
produce the minimum excitation in the 1/akp > (BEC) regime. Thus, this is consistent with the interpretation that 
the studied regime (—1.3 < 1/akp < —0.3) is that of strong three particle correlation. The energies for the broken 
pair states are shown in the Fig. [3] for the akp = —1.3 case. The gap A is estimated from the data sets represented 
by circles with total N = 18, 19, 20 and N = 21, 22, 23 respectively. The calculated A/E FG « 0.3(3). The error bars 
are large for the pairing gap since the three particle effects predominates rather than two particle pairing. In the 
Fig. [3J we notice that the energy dips when total N is a multiple of 3. We interpret this as an effect analogous to 
that observed in the <?(r)'s. This is the evidence that the trimer interaction that brings the green-red pair and the 
blue particle together is much stronger than simple pairwise interaction. Thus, completing green-red-blue trimer is 
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energetically more favorable than unbalanced excess of one or two species. In fact, the trimer binding energy is so 
strong that A < \E trimer /3\ in contrast to A \E pair /2\ of s — 2 Fermi gas in the BEC regime. Consequently, in the 
quasiparticle spectrum, it is expected that we can observe two distinctive gaps; one due to the superfluid pairing and 
another due to the trimer binding. We also notice that at N — 21 the dip is as deep as at N — 24 which indicates 
possible shell closure effect. As seen in the ground state energy (Fig. [1]) and the overlap of the wave functions (Fig. 
|4]), the broken symmetry superfluidity is suppressed in the ak,F > —0.5 regime where the SU(3) symmetry is restored. 

For the comparison purpose only, we naively consider extension of the mean field (BCS-Leggett) method to the 
three component Fermi gas. For this mean field model, there is no R dependence and the variational ground state is 
always stable. The two color pairing gap of the degeneracy three superfluid is given by the same relation as that of 
the degeneracy two Fermi gas. While in the original BCS formalism, the chemical potential fi is kept constant 
and Ascs/ e F = \e*l 2akF we consider the BCS-Leggett [26[ variational formalism {A-BCS-Leggett)- Here, the 
chemical potential is changed in order to keep the density constant. The condensation energy of the system is 
E CO nd — — n o%" where n is the state density (n = Y^fM ) . Dependence on the degeneracy s is included in the 
condensation energy per particle Ec ™ d = —j^-j^: because of the relation ^f-p = hp. Thus, at l/akp — —1, we 

estimate that ^BCS-Leggett/ Efg = 0.33 and E con d/(NEFG) ~ —0.016. At \jakp = —0.75 (or akp = 1.3) ,we have 
^BCS-Leggett/EpG = 0.5 and E con d / {N E p g) ~ —0.038. This estimate is close to the one calculated by GFMC (Fig. 
[3l A/Epc ~ 0.3 ) at the same interaction strength. 

Although for s = 3 Fermi gas, R dependence cannot be removed, clear qualitative differences between the s = 2 
and s — 3 Fermi gases emerge. Unlike in the s — 2 Fermi gas where both the superfluid pair and the bound state 
are qualitatively similar, in the s = 3 Fermi gas the paired states decouple from the bound state(trimer) in energy. 
This can be clearly observed in the quasiparticle excitation spectrum. Realistic interaction potential and channel 
dependence of the interactions are necessary in order to produce not only qualitative but also quantitatively correct 
results for a given three component Fermi gas. We found a regime of interaction strength where the broken symmetry 
pairing is clearly detectable, beyond which three particle effects dominate and the symmetry is restored. Also, non 
trivial dependence of the <?(?") 's on the pairing correlations was discussed. This work has been supported in part by 
the US National Science Foundation via grant PHY 00-98353 and PHY 03-55014. SYC also acknowledges support 
by the DARPA grant BAA 06-19. The authors acknowledge useful comments from J. Carlson and A. Bulgac. One 
of the authors(VRP) passed away during the preparation of the present manuscript and the work was posthumously 
completed. 
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